Moving Wigner Glasses and Smectics: Dynamics of Disordered Wigner Crystals 
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We examine the dynamics of driven classical Wigner solids interacting with quenched disorder 
from charged impurities. For strong disorder, the initial motion is plastic — in the form of crossing 
winding channels. For increasing drive, the disordered Wigner glass can reorder to a moving Wigner 
smectic-with the electrons moving in non-crossing ID channels. These different dynamic phases can 
be related to the conduction noise and I(V) curves. For strong disorder, we show criticality in the 
voltage onset just above depinning. We also obtain the dynamic phase diagram for driven Wigner 
solids and prove that there is a finite threshold for transverse sliding, recently found experimentally. 



PACS numbers: 73.50.-h 

Ordered arrays of charged particles have been studied 
in the context of colloidal suspensions, ion rings, atomic- 
ion Wigner crystals, quantum computers, astrophysics, 
biophysics, plasmas, electrons deposited on liquid He- 
lium surfaces, semiconductor heterostructures, Wigner 
droplets, and arrays of metallic islands interconnected 
by tunnel junctions (!]]. A revival of general interest in 
charged arrays has been fueled by the observation, in 
2D heterostructures, of nonlinear I-V curves exhibiting a 
threshold as a function of an externally applied electric 
field |§JjJ, indicating the presence of a Wigner solid (WS) 
H that has been pinned by the disorder in the sample. 
The depinning threshold can vary up to two orders of 
magnitude in different samples ^j. fndeed, experiments 
jjj based on transport and photoluminescence provide 
indirect evidence (without observing the location of the 
charges) for the existence of the WS, and demonstrate the 
very important role disorder plays in the dynamics of the 
WS || . It is the purpose of this work to study how dis- 
order affects the transport properties of the driven WS. 

The onset of broadband conduction noise has been in- 
terpreted as a signature of the sliding of a defected WS 
If the electrons were to retain their order and slide 
collectively, narrow band noise resembling that seen in 
sliding charge-density wave systems should be observed 
instead. Simulations by Cha and Fertig on classical 
Wigner crystals interacting with charged defects indicate 
that a disorder-induced transition, from a clean to a de- 
fected WS, can occur as a function of increasing pin- 
ning strength jg]. For strong pinning, the initial depin- 
ning is plastic and involves tearing of the electron crystal 
H (b) . However, many aspects of this transport have not 
been systematically characterized, including the current- 
voltage I(V) curves, conduction noise, transverse mean- 
dering, transverse threshold for sliding, and the electron 
lattice structure and motion for varying applied drives. 

Similar plastic depinning transitions have been ob- 
served in the related system of driven vortex lattices in 
disordered superconductors. The plastic depinning of the 
vortices is associated with flux motion through intricate 
river-like channels |R,p3. Defects in the vortex lattice 



strongly affect the depinning thresholds and the voltage 
noise signatures produced by the system. For increasing 
drive the initially defected vortex lattice can reorder to 
a moving lattice or a moving smectic phase p-|Io||. In 
the fast moving lattice phase the vortex lattice regains 
order in both transverse and longitudinal directions with 
respect to the driving force, while in the smectic state 
only order transverse to the driving direction appears. 

It is unclear a priori whether the same type of reorder- 
ing transitions can occur in a system containing long- 
range pinning, such as in the case of the WS interacting 
with charged impurities. This is in contrast with the 
vortex system, where the pinning interacts with the flux 
lines only on a very short length scale. Additionally, 
in the strong pinning limit, critical behavior may occur 
near the depinning threshold, leading to velocity-force re- 
lations |l3| of the form v ~ (F — Ft)^- For instance, for 
transport in metallic dots, £ = 5/3 theoretically |Q, and 
£ = 2.0 and 1.58 experimentally p"5[ . 

In this work, we use numerical simulations to study 
dynamical transitions in a 2D electron system forming a 
classical WS in the presence of charged impurities. For 
strong pinning, the electron crystal is highly defected and 
depins plastically, with certain electrons flowing in well- 
defined channels or interconnecting rivers, while others 
remain immobile. For increasing applied driving force, 
the initially disordered electrons can partially reorder to 
a moving Wigner smectic state where the electrons flow in 
ID non-crossing winding channels. The reordering tran- 
sition is accompanied by a saturation in the dl/dV curves 
and by a change from a broad-band to a narrow-band 
voltage noise signature. For weak disorder, the depinning 
is elastic and a narrow-band noise signal can be observed 
at all drives above depinning. For strong disorder, where 
the depinning is plastic, we find criticality in the velocity 
force curves in agreement with transport in metallic dots 
@,^5| ■ We map out the dynamic phase diagram as a 
function of disorder strength and applied driving force. 
Also, we find unequivocal evidence for a transverse finite 
threshold for sliding conduction, in agreement with re- 
cent experiments 17 . Our results can also be tested in 
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other systems listed in the first sentence of this paper 
(e.g., with charged colloidal systems suspended in water 
in the presence of a disordered substrate). 

We conduct overdamped molecular dynamics (MD) 
simulations using a model similar to that studied by Cha 
and Fertig 0. The energy from interactions is 
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The first term is the electron-electron (Coulomb) repul- 
sion and the second term is the electron-impurity interac- 
tion, where the impurities are positively charged defects 
out of plane. Here, is the location of electron i, and 
is the in-plane location of a positive impurity which 
is located at an out-of-plane distance d (measured here 
in units of do, the average lattice constant of the WS). 
The number of electrons N equals the number of impu- 
rities N p , and the disorder strength is varied by chang- 
ing d. We have also considered many other cases (e.g., 
± charged impurities, Ni ^ N p , etc.) with consistent 
results. The long-range Coulomb interactions are evalu- 
ated with a fast converging sum technique by Jensen ]T^ ] 
which is computationally more efficient than the Ewald- 
sum technique used in previous studies Q . 

To obtain the initial electron positions, we perform 
simulated annealing in which we start from a high tem- 
perature above the electron lattice melting transition and 
slowly cool to a low temperature. Temperature is imple- 
mented via Langevin random kicks. Once the electron 
configuration has been initialized, the critical depinning 
force is determined by applying a very slowly increasing 
uniform driving force fd which would correspond to an 
applied electric field. After each drive increment, we wait 
10 4 MD time steps before taking data, which we average 
over the next 10 4 time steps. For each drive increment, 
we measure the average electron velocity (oc current) in 
the direction of drive, V x = (1/iVi) J2i " x -v i- The V x ver- 
sus fd curve corresponds to an I(V) experimental curve 
and we will thus use the notation V x = I and fd = V . 
The depinning force f2 is defined as the drive fd at which 
V x reaches a value of 0.01 that of the ohmic response. We 
have studied system sizes Ni from 64 to 800 and find sim- 
ilar behavior at all sizes. Most of the results presented 
here are for systems with Ni = 256. 

When driven through a sample containing strong pin- 
ning, the electron lattice undergoes a gradual reorder- 
ing transition as the driving force is increased. We illus- 
trate such a reordering transition in Fig. by plotting the 
eastbound electron trajectories at different driving forces 
for a sample containing strong disorder of d = 0.65. In 
Fig. 0(a) the onset of motion occurs through the open- 
ing of a single winding channel. Electrons outside of the 
channel remain pinned, and the overall electron lattice is 
disordered. At fd/ fd = 1-5, shown in Fig. 0(b), several 
channels have opened, some of which are interconnecting. 



(a) 








(e) (f) 
FIG. 1. Electron positions (dots) and eastbound trajecto- 
ries (lines), for a sample with Ni/N p — 1.0 and d — 0.65. 
fd/fd = ( a ) 1-1. ( b ) 1-5, (c) 2.25, (d) 3.0, (e) 4.0, (f) 5.0. 

The original channel in Fig. 0(a) has grown in width, but 
regions of pinned electrons are still present. Electrons 
moving past a pinned electron perturb it, causing it to 
move (like a revolving door) in a circular orbit around 
the center of the potential minima in which it is trapped. 
Several of these electron turnstiles can be observed in 
Fig. 0(b). In Fig. |(c) for f d /f d - 2.25 there appear 
to be regions where electrons do not flow, but none of 
the electrons are permanently pinned. Some electrons 
become pinned for a period of time before moving again. 
If the trajectories are drawn for a sufficiently long time, 
the electron flow appears everywhere in the sample, al- 
though there are still preferred paths in which more elec- 
trons flow. In Fig. 0(d) for fd/fd = 3-0 the electron flow 
occurs more uniformly across the sample. In Fig. 0(e) 
a t fd/fd — 4.0 the electrons begin to flow predomi- 
nantly in certain non-crossing channels, although some 
electrons jump from channel to channel. In Fig. 0(f) for 
fd/fd = 5.0 the electron flow occurs in well defined non- 
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FIG. 3. (a) I(V) (dashed line) and dl/dV curves (solid 
line) for a sample with d = 0.65. (b) Fraction P$ of sixfold 
coordinated electrons, as a function of driving force fd- (c) 
Fraction Dtr of transversely wandering electrons. Inset to 
(c): Noise spectra for fd/ fd = 1-5 in the plastic flow regime 
showing broad-band noise (dashed line) and for fd/ fd — 4.0 in 
the smectic regime showing a narrow-band signal, (d) Log-log 
plot of V x = {fd - fdf, in which £ = 1.61 ± 0.10 and 1.71. 
(e) V x versus fd for disorder strengths of d = 0.5 and 0.65. 

the degree of order and the nature of the flow in the sys- 
tem. To quantify the degree of plasticity of the electron 
flow, we plot in Fig. [|(c) the fraction of electrons D tr 
that wandered a distance of more than oq/2 in the di- 
rection transverse to the drive during an interval of 8000 
MD steps. The peak in D tr coincides with the peak in 
dl/dV. D tr then slowly declines until it saturates at 
fd/ f% ~ 6, indicating the gradual formation of the non- 
crossing channels as seen in Fig. |l|(f). The saturation in 
D tr coincides with the saturations in both Pq and dl /dV. 

The I{V) curve in Fig. ||(a) corresponds to a highly 
irregular voltage signal as a function of time when the 
electrons are in the plastic flow regime (fd/fc = 1-5). 
The corresponding voltage noise spectrum in the inset of 
Fig. ||(c) shows that only broad-band noise is present. In 
contrast, at fd/ fd = 4.0, in the moving smectic regime, 
a roughly regular signal is obtained, and narrow-band 
noise appears, as shown in the inset of Fig. ||(c). In a 
system with d = 1.57 when the depinning is elastic only, 
an even more pronounced narrow band noise signal is 
observable. In experiments, broad band noise has been 
observed above depinning [Q] , but narrow band noise has 
not been seen. Our results suggest that plasticity may be 
playing an important role in most experiments. 



FIG. 2. Delaunay triangulation for electrons in a sample 
with (a) d — 0.65 (strong pinning) at fd/ fd — 5-0; and (b) 
d — 1.76 (weak pinning) and fd/ fa — 2-0. Large circles indi- 
cate (5- or 7-fold) defects in the electron lattice. 

crossing channels, which can contain different numbers 
of electrons. A similar channel motion exists for driven 
vortices in disordered superconductors pL|l^] . 

For samples containing very weak disorder, the pinned 
WS has six-fold ordering and depins elastically, without 
generating defects in the lattice. In this regime, the 
electron crystal flows in ID channels with each chan- 
nel containing the same number of electrons. Here, the 
transverse wanderings of the electrons are considerably 
reduced compared to the case of strong pinning. 

To better illustrate the change in the amount of disor- 
der in the electron lattice, in Fig. ||(a) we show the Delau- 
nay triangulation for the electrons for the same drive as 
in Fig. |l](f). Defects, in the form of 5-7 disclination pairs, 
appear with their Burgers vectors oriented perpendicular 
to the direction of the drive. In the structure factor there 
are only two prominent peaks for order in the direction 
transverse to the drive, consistent with a moving Wigner 
smectic state. In Fig. |](b) we show the Delaunay trian- 
gulation for the moving state in a weakly pinned sample 
with d = 1.76 where the initial depinning is elastic. Here 
the moving lattice is defect free. The structure factor in 
this case shows four longitudinal peaks in addition to the 
two transverse peaks, although the transverse peaks are 
more prominent. Much larger systems would be neces- 
sary to determine whether the system is in a smectic state 
or in a moving Bragg glass in this case of weak disorder. 

In order to connect the reordering sequence with a 
measure that is readily accessible experimentally, we plot 
in Fig. ||(a) the I(V) and dl/dV curves. A peak occurs 
in dl/dV at fd/ fd ~ 3, when the electrons are undergo- 
ing very disordered plastic flow. In Fig. ||(b) we plot the 
fraction Pq of six-fold coordinated electrons as a function 
of drive. A perfect triangular lattice would have Pg = 1. 
For drives fd/ fd < 2 the lattice is highly defected. For 
fd/ fd > 3 (i.e., past the peak in dl/dV), the order in 
the lattice begins to increase. The value of Pq saturates 
near fd/ fd ~ 6 which also coincides with the saturation 
of the dl/dV curve. Thus the experimentally observable 
I(V) characteristic can be considered a good measure of 
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FIG. 4. Dynamic phase diagram for the driven disordered 
Wigner solid. Inset: Clear evidence for a finite transverse 
depinning threshold for a system with d = 0.9 at fd = 0.16 in 
the reordered phase. 

In Fig. ||(d,e) we examine the critical behavior above 
the depinning threshold of V x versus fd for disorder 
strengths d = 0.5 and 0.65. In Fig. ||(d), a log-log plot 
of V x ~ [(fd — f%)/fffi, indicates that the curves are 
fit well by a power law over one and half decades with 
£ = 1.61 ± 0.10 and 1.71 ± 0.10, respectively. These val- 
ues agree well with the predicted value of £ = 5/3 
for electron flow through disordered arrays, and the ex- 
perimentally observed values of £ = 1.58 and 2.0 Jlj|. 

The theory of the moving ordered phase || predicts 
a barrier to transverse motion once channels similar to 
those in Fig. 1(e) form. As shown in the inset of Fig. 4, 
for a system in the reordered phase we observe a trans- 
verse depinning threshold that is about 1/6 the size of 
the longitudinal depinning threshold /J. In recent ex- 
periments, Perruchot et al. jl7| also find evidence for a 
transverse barrier that is about 1/10 the size of the longi- 
tudinal threshold. These thresholds are much larger than 
those observed in vortex matter interacting with short- 
ranged pinning, where ratios of 1/100 are seen [ jll| . 

In Fig. [| we present the WS dynamic phase diagram 
as a function of disorder strength and driving force. For 
d < 1.4 there are a considerable number of defects in the 
WS, and the initial depinning is plastic. We label this 
region the pinned Wigner glass. For increasing disorder, 
the pinned region grows while for strong enough drive 
the electrons reorder to a moving Wigner smectic state. 
For d > 1.4 there are few or no defects in the pinned 
state and the initial depinning is elastic. The presence 
of a crossover from elastic to plastic depinning with de- 
creasing disorder strength is in agreement with §. 

In summary, we have investigated the pinning and dy- 
namics of an electron solid interacting with charged dis- 
order. We find that for strong disorder the depinning 



transition is plastic with electrons flowing in a network 
of winding channels. For increasing drives, the electrons 
partially reorder and flow in non-crossing channels form- 
ing a moving Wigner smectic. We show that the onset of 
these different phases can be inferred from the transport 
characteristics. In the plastic flow regime, the noise has 
broad-band characteristics, while in the moving smec- 
tic or elastic flow phase, a narrow band noise signal is 
observable. We also show that the onset of the plastic 
flow phase shows critical behavior with critical exponents 
in agreement with predictions for transport in arrays of 
metallic dots |l4|,[l^]. We map out the dynamic phase 
diagram as a function of disorder and applied driving 
force. We obtain a finite threshold for transverse sliding, 
in agreement with recent experiments (ItJ . 
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